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The radiative transfer equations for multiple inverse Compton scattering of the Cosmic 
Microwave Background Radiation (CMBR) by the hot intra-cluster electrons are solved 
numerically. The spherical isothermal and inhomogeneous j3 model has been considered 
for the electron distribution. The anisotropy of the CMBR caused by scattering, known 
as thermal Sunyaev-Zel'dovich effect, along the radial axis of the medium is compared 
with the analytical solution of Kompaneets equation. The X-ray data of several clusters 
of galaxies at low redshifts provide an estimation of the central electron density no to be 
of the order 10 -3 . It is found that for this value of no the effect of multiple scattering is 
negligible. The numerically calculated anisotropy along the radial axis matches well with 
the analytical solution that describes single scattering. The result incorporating multiple 
scattering is fitted with the recent observation of Sunyaev-Zel'dovich effect in the cluster 
Abell 2163. It is shown that if no is greater by an order of magnitude, which could be 
possible for cluster of galaxies at comparatively higher redshift, multiple scattering would 
play a significant role at the Wien region of the anisotropy spectrum. A fitting formula 
for the correction to the Sunyaev-Zel'dovich effect due to multiple scattering is provided. 
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1. Introduction 

Inverse Compton scattering of the Cosmic Microwave Background Radiation 
(CMBR) by hot intra-cluster gas - better known as the thermal Sunyaev-Zel'dovich 
effect (TSZE) 1-3 results in a systematic transfer of photons from the Rayleigh-Jeans 
to the Wien side of the spectrum causing a distortion in the Planckian nature of 
the spectrum. The measurements of the effect yield directly the properties of the 
hot intra-cluster gas, the total dynamical mass of the cluster as well as the indirect 
information on the cosmological evolution of the clusters. The effect is also used as 
an important tool to determine the Hubble constant Hq and the density parameter 
fio of the Universe 4 ^ 10 . Recent interferometric imaging of the TSZE n ~ 14 has been 
used to estimate the mass of various galaxy clusters which could constrain the cos- 
mological parameters of structure formation models. Measurements of TSZE and 
the kinetic Sunyaev Zeldovich effect (KSZE) from the Sunyaev-Zel'dovich Infrared 
Experiment and OVRO/BIMA determine the central Compton y parameter and 
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constrain the radial peculiar velocity of the clusters. In near future sensitive ob- 
servations of the effect with ground based and balloon-borne telescopes, equipped 
with bolometric multi-frequency arrays, are expected to yield high quality measure- 
ments. 

In modeling the observational data, the analytical solution of Kompaneets 
equation 15 with relativistic corrections 16 are used in general. All theoretical discus- 
sions and observational inferences so far are restricted to the case of single Compton 
scattering. Single scattering approximation is well justified because X-ray data of 
most of the galaxy clusters wherein the Sunyaev-Zel'dovich effect has been inves- 
tigated, indicate low optical depth of the intra-cluster medium (ICM). Effects of 
multiple scattering have been considered by many authors 17-21 . Recently, Itoh, 
Kawana, Nozawa and Kohyama 22 have studied the effect of multiple scattering by 
performing direct numerical integration of the collision term of the Boltzmann equa- 
tion and found that the effect of multiple scattering is negligible for the observed 
galaxy clusters. 

In this paper, the numerical solutions of the most general radiative transfer 
equation for inverse Compton scattering of CMBR photon are presented by incor- 
porating multiple scattering in a spherically symmetric inhomogeneous and isother- 
mal medium. Although my approach is different, the results are in good agreement 
with that of Itoh et al. 22 . I further present the results by considering a medium 
with optical depth an order of magnitude higher than that of the galaxy clusters 
observed so far. 

In section 2, 1 present the most general radiative transfer equation that describes 
the inverse Compton multiple scattering of low energy photons with non-relativistic 
electrons in an inhomogeneous medium. The Kompaneets equation which is a spe- 
cial case of isotropic scattering is also presented with the analytical solution that 
describes the anisotropy in the CMBR spectrum due to single scattering. The model 
parameters are provided in section 3. The results are discussed in section 4 followed 
by the conclusion in section 5. 

2. The Radiative Transfer Equations 
2.1. General case 

Chandrasckhar 23 has provided the radiative transfer equations for the Compton 
scattering by assuming that the electron energy is much less than the photon energy. 
On the other hand, for the scattering of CMBR, the photon energy is much less 
than the electron energy. In fact, in the TSZE, the electrons are considered to have 
relativistic motion described by relativistic Maxwellian distribution. The relevant 
radiative transfer equation in its general form that describes the inverse Compton 
multiple scattering of low energy photons in a spherically symmetric inhomogeneous 
medium can be written as 24-26 




dl v 



dl v {fi,r) l-^ ; 



dr r 



? dl v (n,r) 

9yU 
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= -a T n e (r)I v (fi,r) + ^<r T n e (r)uj a x 

f U - P + 2upi P 2 )lAp,rW. (1) 

J_ 1 [ m e c z m e c z \ ov z av ) J 

The above equation describes the inverse Compton scattering of low-frequency 
radiation (first order in hv/mc 2 ) in a hot thermal isotropic electron gas (first order 
in kT e /mc 2 ) in the rest frame of the electron gas. Here loq is the albedo for single 
scattering, s is the ray path, /j = cos 8 where 9 is the angle between the axis of 
symmetry (radial axis) in the rest frame of electron gas and the ray path, ar and 
n e being the Thomson scattering cross section and the electron number density 
respectively, fc, c, m ei v and T e are Boltzmann constant, velocity of light, electron 
rest mass, frequency of the photon and the temperature of the electron respectively. 
coo = 1 for a purely scattering medium. The phase functions Pq, Pi and Pi are given 
as : 

P (m,m') = ^[3-M 2 -m' 2 (1-3m 2 )], (2) 
Pi{n, (J,') = | [1 - 3 M ' 2 - 3^ 2 (1 - 3/i' 2 ) + 2 M V (3 - 5m' 2 ) + 2mm'(3m' 2 - 1)] , (3) 



and 

P 2 (M, m') = | [3 - M 2 - M' 2 + MM' (3MM' - 5 + 3m 2 + 3m' 2 - 5mV 2 )] • (4) 



2.2. Special case : Isotropic radiation field 

For an isotropic radiation field 

\ j /(/i'.l/Jd/i' =/(!/), 

y j(m»m'<v = y i(m,^)m' 3 ^' = o. 

Therefore, for isotropic scattering, equation (1) reduces to 



and 



(5) 
(6) 

(7) 



— =a T n e (r)— 5 i/ -5-=- -2^— (8) 



which is well known as the Kompancets equation for low frequency photons. 

Integration of equation (8) along the ray path provides the thermal component 
of the distortion AI and is written as 

x^e x r x(e x + 1) 



(1 + Jt) (9) 
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where x = kT ^ BR , Io = 2(kT CM BR) 3 /(he) 2 . The term y = ^ J <j T n e ds is 



usually referred to as the Compton y parameter and St is a relativistic correction 
to the thermal effect 16 significant if kT e > 5KeV. 

3. The Models 

I have adopted the isothermal model 27 for the density distribution of the clusters. 
The spherical model density is described by 



where the core radius r c and are shape parameters, no is the central electron 
density. 

I have adopted the values of the model parameters given in LaRoque et al. 13 for 
the Sunyaev-Zel'dovich anisotropy spectrum observed in Abell 2163. A combination 
of ROSAT X-ray data and OVRO/BIMA observation of Sunyaev-Zel'dovich effect 
of Abell 2163 at the cosmological rcdshift z=0.203 implies the electron temperature 
T e = 12.4 KeV, the angular core radius of the cluster 9 C = 1'.20 ± O'.ll and the 
shape parameter = 0.616±0.031. However, for simplicity, I have adopted = 2/3 
which also provides a good fit to the observed data. Following LaRoque et al. 13 , I 
have considered a cosmological model with the matter density parameter Q m = 0.3 
and the cosmological constant corresponding to fl v = 0.7. I have taken the Hubble 
constant Ho = 58.0 Kms _1 Mpc _1 . With the above parameters, observational fit 
of the Sunyaev-Zel'dovich effect for Abell 2163 yields the central electron number 
density no — 6.96 x 10~ 3 and hence the corresponding Compton y parameter is 
3.83 x 10~ 4 . In order to investigate the effect of multiple scattering in a denser 
medium I have also considered models with 10 x no and 20 x no keeping other 
parameters unaltered. 

4. Results and Discussion 

Equation (1) is a coupled integro partial differential equation and so it cannot be 
solved analytically. I have solved it numerically by discretization method. In this 
method the medium is divided into several shells and the integration is performed 
over two dimensional grids of angular and radial points. The numerical method 
is described in detail in Peraiah & Varghese 28 . I have used 11 points Trapezoidal 
method for angular and frequency integration. The numerical code is thoroughly 
tested for stability and flux conservation. 

For the boundary condition, I have provided equal amount of intensity corre- 
sponding to Tcmbr = 2.728K at optical depth r = along all directions, i.e., 




(10) 



J(±/i, v, r = 0) = (eWWbMB* _ i) 



(11) 



The relativistic correction (1 + <jt) is included to the emergent intensity. 
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Fig. 1. Best model fit for thermal Sunyaev-Zel'dovich spectrum of Abcll 2163. Solid line repre- 
sents the distortion along the radial axis for multiple scattering and dashed line represents that 
along the ray path for single scattering. The central electron number density no = 6.96 X 10 -3 



Equation (1) decribes the anisotropic and inhomogencous radiation field 
in a spherically symmetric medium. Although scattering makes the radiation 
anisotropic, spherical symmetry of the medium and the initial isotropy of the radi- 
ation make the emergent radiation field isotropic. As a consequence the distortion 
along the radial axis and that along the ray path is the same. 

The numerical results that incorporate multiple scattering and the results for 
single scattering are presented in Fig 1, Fig 2 and Fig 3 for different values of the 
electron number density. Since the optical depth with a central electron number 
density n = 6.96 x 10~ 3 is small, the effect of multiple scattering is negligible 
as seen in Fig 1. The slight differences in the results at the Wien and Raylcicgh 
Jeans region is well in agreement with that obtained by Itoh et al. 22 . This can 
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Fig. 2. Same as Fig. 1 but with the central electron density 10 X no- 

be visualized from Fig. 1 where I have fitted the observational data for the galaxy 
cluster Abell 2163. The model fit clearly shows that the effect of multiple scattering 
can be neglected for Abell 2163 within the uncertainties in the physical parameters. 

In Fig. 2 and Fig. 3 the effect of multiple scattering is shown by increasing the 
central electron density ten and twenty times respectively. The other parameters 
such as the electron temperature, shape parameters etc. and the cosmological pa- 
rameters are kept unchanged. With the increase in the electron number density 
and hence the optical depth of the medium, the number of scattering increases. 
Consequently, the distortion in the spectrum increases significantly. The increase 
in the distortion is much more at the Wien region of the spectrum than that at the 
Rayleigh- Jeans region. 

The spectral distortion due to Sunyaev-Zel'dovich effect by single scattering is 
characterized by three distinct frequencies : the cross over frequency xq = 3.83 
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Fig. 3. Same as Fig. 1 but with the central electron density 20 X no- 
where the thermal Sunyaev Zel'dovich effect vanishes, x m i n — 2.26 which gives 
the minimum decrement of the CMB intensity and x max — 6.51 which gives the 
maximum distortion. The value of x n is pushed to higher values of x with the 
increase in T e for the relativistic case. Fig. 2 and Fig. 3 shows that there is no chnage 
in the values of xq , x m i n and a negligible change in x max when multiple scattering is 
included. The amount of distortion due to multiple scattering is however higher at 
Xmin as compared to that with single scattering. At x max the amount of distortion 
remains almost the same even if the central electron density is increased by twenty 
times. The effect of multiple scattering becomes significant as x increases from x max 
that is at the far Wien region. Therefore, multiple scattering must be incorporated 
in modeling the thermal Sunyaev-Zel'dovich effect anisotropy spectrum if the optical 
depth of the medium is comparatively high. Quantitatively, if the central electron 
density of a galaxy cluster is one order more than that of Abcll 2163 or other clusters 
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observed so far then the effect of multiple scattering would be very important. 

Finally, the thermal component of the distortion AI modified due to the multiple 
scattering can be written as : 

~x(e x + 1) 



AI(x) = I Q y 



x 4 e x 



(e x - l) 2 



(l + *r)(l + 5 



AI J 



(12) 



e x - 1 

where the correction due to multiple scattering 5m is expressed by a fitting formula 
given by 

'0.139 s " 



7.2n 



1 + log 



n 

-2„,4 



(-1.864897 + 2.744115x - 1.337648x 2 + 0.307672a; J 



-5^6 



-3.812115 x 10^x 4 + 2.634191 x HTV - 9.551485 x lO'V 
1.432487 x 10" V). 



(13) 



5. Conclusions 

The general transfer equations for inhomogeneous low frequency radiation multiply 
scattered by relativistic thermal electron gas in a spherically symmetric; medium are 
solved numerically. The observed Sunyaev-Zeldovich effect for Abell 2163 is fitted 
with the numerical result. The numerical result for the distortion in the CMBR 
due to multiple scattering is compared with that given by the analytical solution of 
Kompaneets equation. The effect of multiple scattering is found to be insignificant 
for the inferred optical depth of the galaxy clusters observed so far. However, if 
the central electron number density is increased by an order of magnitude, multiple 
scattering affects the CMBR distortion spectrum significantly in the Wien region. 
Therefore, it is important to include multiple scattering while modeling the observed 
Sunyaev-Zel'dovich effect of galaxy clusters with high electron density. A fitting 
formula for incorporating the correction due to multiple scattering is provided. 
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